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ABSTRACT 

A reexamination of Hugoniot equation of state data and three new release 
adlabat points indicates that results for enstat lte-bronzite corposltion 
pyroxene are compatible with its transforming to a perovskite phase at high 
pressure (for En^o 5 PO ■ 4.20 gm/cm\ K 0 ^260 ± 35 GPa, ^3.5 i .65). 

The release adlabat data as well as results from porous samples Imply 
that the shock-wave data do not define an equilibrium, high pressure 
phase Hugoniot below about 100 GPa. These also suggest a further 
transformation to a phase (or assemblage) with density about 5% (or more) 
greater than that of orthorhombic perovskite. The data would allow such 
a transformation to occur at pressures as low as '^60 GPa under shock, 
representing an upper bound for the equilibrium transition pressure. 

Hugoniot data on magnesian olivines also appear to represent states ol 
thermodynamic disequilibrium or a mixed phase region below about 80-100 GPa. 
However, Hugoniot points for Mg-pyroxene and Mg-olivine coincide at pressures 
above 70 GPa suggesting that these minerals transform to high pressure phases 
or phase assemblages of comparable density. Since MgO (presumably as periclase 
or in a closely related structure) attains relatively low densities at these 
pressures, the shock-wave data are in strong disagreement with the dispropor- 
tionation of Mg2Si04 to a MgSi03 (perovskite) + MgO assemblage above 80-100 GPa. 
Conversely, the shock-wave data do not preclude a transformation of the type, 
Mg2Si20g -+ Mg2S J .04 ("post-perovskite" phase) + Si02 (rutile or fluorite struc- 
ture). Again, these results would Imply polymorphism to very dense, "p«st- 
perovskite" phases. Based on the arguments made for pyroxene, this ould 
occur at pressures as low as 60 GPa under equilibrium conditions. 

We note that the combined results of high-pressure experiments allow 
Mg-pyroxene compositions to be as likely candidates for the lower mantle as 


olivine. 
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INTRODUCTION 

Upon recognition that the features associated with plate tectonic processes 
on the earth (subduction zones, ridges and plumes) are either minor surficial 
expressions, or alternatively reflect an integral sampling of the mantle con- 
vective flow field, and with the perspectives attained from the last decade 
of terrestrial planetary exploration, experimental high-pressure geophysics 
has focussed increasingly on studying candidate materials of the earth's lower 
mantle and core. 

Current ideas (e.g., Schube rt and Ande rson , 1975] suggest that active 
convection in the earth takes place to depths of at least ^750 km. However, 
it is important to discern whether the upper mantle above ^350 km depth, and 
possibly the phase transition region extending to depths of at least M)70 km, 
represents a differentiate of the lower mantle and is thus atypical of the 
silicate zone of the earth as proposed by D. L. Anderson et al. [1972] and 

ID. L. Anderson [1976], in contradistinction to the views of Ringwood [1975] 

which have emphasized the essential homogeneity of the mantle. 

Clear-cut data specifying the composition of the upper mantle stems 
largely from inclusions in kimberlites [e.g., Boyd , 1973] and nodules from 
volcanic rocks of deep-seated origin [’lacGrcgo., 1975 and Basu, 1975. for 
example] which demonstrate the upper mantle to be largely garnet lherzolite 
with minor amounts of eclogite. This composition is not unlike the "pyrollte" 

model originally proposed by Clark and Ringwood [1964] which, by means of varia- 
tions in temperature, pressure, and degree of partial melting, accounts for the 
genesis of a large class of effusive and plutonic rocks which penetrate the crust 
[ Green , 1972]. 

The lower mantle, extend* ■g from depths of 670 to 2900 km, represents 
49% of the mass of the earth. What can be learned of its mineralogy and 
thermal state must, at present, be inferred from comparison of shock-wave 
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m<l MlMO logical tl.it. i [Al'tschulyr .uni Sharp ldzhanov, 1971 . 1 , !> ; Amlrrsun 
et al . , 1972; and Andergon 1976], and by extrapolation of results based on 
quench products from static high pressure, petrological experimentation 
{e.g. Rlngwood , 1975; Aklm oto et al . , 1976; Liu, 1975a, b, 1976] using largely 
ultrasonic and static compression equations of state. 

Not unrelated to the questions of the properties of pyroxenes and olivines 
at high pressures and hence, of mantle compositions, has been an increasing 
awareness among earth scientists, that the composition, structure, and thermal 
history of the earth must ultimately be related to the composition of the sun, 
the moon, meteorites and other terrestrial planets, and the processes of 
planetary accretion (e.g. L ewis , 1972; Rl ngwood , 1975; Hanks a n d Ander son, 19691. 

In the present study the nature of the shock- induced, high-pressure phases 
of olivine and pyroxene rocks is examined in the light of data for the 
densities of a new class of perovskite-.elated silicate structures synthesized 
by Liu [1975a, b, and 1976], and some new Hugoniot and release adiabat data for 
bronzite. Previous data for olivine-rich samples in the range FolOO to Fo85 
and for enstatitic samples (over a similar range of Fe contents) are reported 
by Tr unin et al. [1965] and McQueen et al. [1967]. Trunin et al. [1965] and 
recently Sla uikov a nd True in [1973] have reported single datum points for 
four different composttons, extending our knowledge of Hugoniot states for 
ferromagnesian silicates to pressures in excess of 240 GPa. 


EXPERIMENTAL METHODS 


Discs (1 cm-diameter) of Bumie bronzlte (Mg >8tl , PemJSlOj were cored 
from the same aliquot previously studied to shock pressures of 48 GPa ( Ahrens 
and Gaffney , 1971] and machined to thickness tolerances of ''<10.05 mm. The 
density of each disc was determined at least twice by weighing in air and 
In reagent grade toluene (temperature corrections from Be rman, 1939). The 
samples were mounted on 0.5 mm thick Cu or Ta driver plates and Impacted by 
fl/er-plate bearing projectiles at speeds ranging from 5.4 to 6.1 km/sec with 
a light-gas apparatus qualitatively similar to that described by Jones et al. 
[1966], Pure (commercial-grade) Cu and Ta flyer plates, 2.5 mm thick, with 
a minimum diameter of 17 mm (Fig. 1) are hot-press welded Into lexan pro- 
jectile blanks using techniques developed by A. Mitchell (private communica- 
tion, Lawrence Livermore Laboratory). After molding and stress relieving, 
a polyethylene, gas sealing rear portion of the projectile was press fit onto 
the 6 mm long lexan flyer plate (front) portion of the projectile. The final 
projectile diameter was machined so as to provide a ''<0.03 mm Interference 
fit with the light-gas gun launch tube. The resulting projectiles had masses 
of 14 and 20 gm for Cu and Ta flyer plates, respectively. The light-gas gun 
used to accelerate these projectiles employs a 20 m-long, 16 cm diameter 
pump tube in which a ^20 Kg high-density polyethylene piston is used to 
compress the propelling gas, H 2 , initially at a pressure of 0.1 MPa. The 
polyethylene piston achieves speeds of M1.4 km/sec using approximately 3 Kg 
of double-based propellant. The projectile is launched upon bursting of a 
diaphram (at a gas pressure of ''<0.07 GPa) at the breech of the 25 mm diameter 
7 mm long launch tube. The total flight of the projectile from the muzzle 
of the launch tube to the target is ^500 mm and occurs in a vacuum of 10 to 
30 mtorr. Projectile velocities are measured using the flash X-ray method 
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of Jones et al. [1966), except that Interruption of a continuous X-ray beam 
[Long and Mi tche ll, 1972) is used to trigger the first flash X-ray unit 
(^15 ns exposure), and the breakage of a 0.05 mm diameter copper wire inter- 
cepting the outer edge of the projectile iH used to trigger the second 

flash X-ray unit. The center lines of the two flash X-ray units are ''-350 mm 

% 

apart and the time Interval between X-ray flashes is measured to within an 
uncertainty of *0.01 usee. Since triggering in our initial experiments was 
not reliable, we have had to rely on the redundancy provided by measuring 
the electronic delays of both flash X-ray units, and the time interval between 
the second X-ray flash and the closure of the redundant pin swltch(es) 
mounted on the sample assembly (Fig. 1) to infer Impact velocity. The ^3° 
projectile tilt observed in Fig. 1, is highly reproducible and the sample 
assembly was thus pre-oriented to achieve nearly normal impact. 

As in the experiments of Ahr ens et al. [1971 J an electronically triggered, 
image-converter streak camera was used to record shock transit-time through 
the sample. However, in the present experiments the onset of illumination 
on shock arrival at the mirrors (presumably arising from the high temperatures 
behind the shock in the glass) was used to record shock arrivals with the 
streak camera, and the arrival of the shock at the mirror free-surfa . . 
observed to result in the immediate loss of illumination (Fig. 2). We infer 
this effect is caused by the disintegration of the specular mirror surface, 
probably upon incongruent vaporization of the glass, and the accompanying 
sudden adiabatic expansion of the vaporized species. In the present experi- 
ments the streak camera was writing at a speed of 25 mm/usec for a duration 
of ^2 usee on 20 X 22.5 cm (4 X 5 inch) Tri-X film. Developing procedures 
were used to force the film to a sensitivity of ASA 3000. Time calibration 
was obtained using a Pockel-cell modulated Ar-ion laser beam modulated at 20 
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Mhz, thus providing time marks at 50 psec intervals. The writing rate at 
any point on the streak camera trace was measurable to within ±0.25%, and a 
practical time resolution of ^2 nsec was achieved. 

Hugoniot states are determined using the impedance match method Rice et al. 
[1958] and release adiabat data employ the buffer technique described by Ahrens 
et al. [1969]. The equations of state for Ta and Cu employed in the impedance 
match solutions are those given by McQueen et al. [1970]. 


7 . 


EXPERIMENTAL RESULTS FOR GLASS AND BAMLF, BRONZITE 

..ie unexpected sel f-i 1 lumlnat Ion and prompt light decay upon passage of 
strong .shocks through the soda-lime glass mirrors (microscope slide quality, 
plate glass) provided consistent, and we believe reliable measurements of 
shock-wave velocities in this material. The resulting data are given In 
Table 1 and plotted in Figures 3 and 4. The shock wave velocity (U ) - 
particle velocity (U^) relation In Fig. 3 Is concordant with data reported 
by Dreatn and Ada durov [1964] to pressures of 41 GPa for a glass of similar 
zero-pressure density (2.48 gm/cm , ) but of slightly different composition. 

The linear U^-U^ relation indicated in Fig. 3 fits the present four data 
points, and Dremln and Adadurov's data points at 33 and 41 GPa with a co- 
efficient of correlation r ■ 0.997. This fit was used to obtain the release 
adlabat points for the Bamle bronzite (Figure 5, Table 2). 

Qualitatively, the behavior of Na-Ca glass is observed to be similar to 
that first reported by Wackerle [1962] for fused quartz. It appears that 
above ''-IS GPa, the Si +4 ion begins to transform from tetrahedral to octahedral 
coordinations with the 0 , '“ ion, as in the case of the quartz to stishovite 
transition [McQ ueen et al ., 1963]. This, and possible coordination increases 
of the Na + an and Ca* 4 ions, account for the lmcompresslble behavior Indicated 
above 50 GPa. The zero-pressure (P 0 ) density shown in Fig. 4 is calculated 
using our microprobe analyses, standard oxide molar volumes [ Robie and W aldbaum, 
1962], and assuming a stishovite molar volume of 14.014 cm 3 for S10 . 

Although the uncertainties in the measurement of projectile and shock 
velocity can be obtained from scatter of redundant projectile time-interval 
measurements and uncertainties in discerning shock arrival-times on the 
streak records, the calculation of uncertainties in the pressure-density 
plane shown in Figures 4 and 5 are less straightforward. The non-orthoganality 
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of the uncertainty brackets shown In both figures for the Hugunlot data 
arise from the fact t'u.t the particle velocity, and hence pressure and density, 
are not Independently related to uncertainties in projectile and shock velocity 
via the Impedance match solution. The uncertainty arlcing from both the 
density and pressure are calculated from the uncertainty in projectile velocity 
holding the shock velocity fixed (at Its average value) and vice versa. The 
uncertainty in the three release points take Into account only the uncertainty 
in shock velocity through the buffer (1 nm-thlck glass). 
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DISCUSSION 

Shock-wave data for magnesian pyroxenes are summarized In Figure 5, 

Including results for porous (synthetic aggregate) samples, non-porous 
samples, and our new data on Bamle bronzlf (cf. Table 2). At pressures 
below about 70-100 CPa It Is clear that the porous Hugonlot data bear an 
anomalous relation with respect to the non-porous data since the former should 
lie on "hot compression" curves, systematically displaced to lower densities 
(at a given pressure) from the principal Hugonlot presumably defined by 
non-porous data. The energy associated with irreversible compaction of 
porosity represents a thermal pressure as a function of Initial density, via 
Grllnelsen's Y (V,T), which defines this displacement. For this reason, the porous 
data have, in general, been very scantily discussed, considered uninterpretable , 
or outright dismitsed as wrong [ McQueen, MarBh, and Fritz , 1967; McQuee n , 1968; 
Ahrens, A nderson, at .d Rtn gwood , 1969; D avi es and Gaffney , 1973; Slmakov and 
Trunin, 1973). 

The relatively high density attained by the porous materials on shocking can, 
nevertheless, be explained as a kinetic effect [c.f. McQueen et al. . 1967. and 

1970). Assuming that pyroxene transforms to a dense phase or assemblage above a 
given pressure, this transformation, even once Initiated under shock, may 
not have time to go to completion without substantial overdriving. Because 
porous samples become much hotter on shocking, their rate of t ransformat Ion 
Is expected to be significantly enhanced, thus yielding a larger proportion 
of high pressure phase material for a given pressure compared with non-porous 
samples. Assuming this model applies, it Is easy to explain the distribution 
of Hugonlot data in the figure. 

The pressure-release points shown in Figure 5 provide independent support 
for such a model. Since the release path is adiabatic (isentropic to within 
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The approximation of reversibility) while entropy increases monotonical ly up 
a material's Hugonlot, the P-V slope of the release path should be no steeper 
than the Hugonlot slope at a given pressure [see DuvalJ_ and Fowles, 1966, for 
example]. Because the measured release adiabat points lie below the Hugonlot 
data, this indicates that a mixed phase region extends to pressures around 
70-100 GPa with transformation to a high density assemblage complete at pressures 
no lower than these. This high pressure assemblage must have a metastable 
Hugonlot with a P-V slope at least as great as the release paths, and thus 
supports the hypothesis of an extended mixed phase (or disequilibrium) region 
probably associated with the rate of phase transformations under shock. 

At pressures above about 100 GPa, the non-porous samples achieve densities 
thut are generally (and consistently' higher than the porous points. So we 
believe that the shock-wave data begin to define an equilibrium, high-pressure 
assemblage Hugonlot at these pressures, below which one or more phase t ransforma- 
tlons result in a wide mixed-phase region. 

Similarly, Hugonlot data for magnesian olivines are summarized in Figure 6. 
These show the same anomalous relations between porous and non-porous samples, 
implying transformation to a high density phase or assemblage subject to 

kinetic effects (though release adiabat data for such compositions are not 
yet available). This conclusion fits in with the discussion of Ahrens a nd 
Peters en [ 1969 ] . 

Aside from these results with porous data and steep measured release 
paths [c.f. Ahrens et a l., 1969; Grady , 1975 and this volume, as well], 
a comparison of shock-recovery and static high-pressure experiments [ Ahrens 
an d Graham , 1972; Liu , 1974 , 1975c; Schne ider and Hor nemann, 1976] also 
Indicate that overdriving beyond equilibrium transition pressures is commonly 
needed in silicates to achieve a given phase transformation under shock, 
clearly implying control by kinetics of nucleation and growth of high pressure 
phases [but note Podurets and Trunin , 1974], Thus, porous samples may well 



11 . 


Indicate the presence and nature of hint) pressure phases (or assemblages) 
more readily than non-porour samples under shock. Also, high pressure phase 
Hugonlots will tend to be significantly steeper (and therefore centered at 
higher Pq densities) than would be apparent from the distribution of the data. 

Zero-pressure densities for high-pressure phases Interpreted from pyroxene 
and olivine shock-wave data by several workers are summarized In Figure 7. 

The tendancy has been to assign higher densities (for a given composition) to 
olivine (A 2 BO 4 ) than to pyroxene (ABO 3 ) stoichiometries. One problem is 
apparent near the magnesian end, however, In that Interpreted densities seem 
to actually decrease with Increasing Fe-content. Barring very complex phase 
relations, this 1l highly unlikely due to the enormous difference in the atomic 
masses of Fe and Mg. Rather, It appears that the density estimates for the 
magnesian endmembers have been biased upwards relative to the Fe-bearlng 
compost tons by the high densities Implied by the porous data (which only 
exist for the endmembers, forsterite and enstatite). Rased on the arguments 
developed above, it is not appropriate to directly eompire the porous data 
with the data for (non-porous) Fe-bearing compositions. 

hiu'a recent (1975a, b, 1976; Liu and Rlng wood. 19751 static high pressure 
experiments have confirmed the previously suspected occurrence and significance 
of perovskite-related structures in silicates [ Reid and Rlngw ood, 1975, 1974, 
1969; Rlngwood and Major , 1971; Shim izu et a l . , 1970). Distorted (non-cublc), 
perovBkite- 1 ike structures have been quenched from diamond-anvil experiments 
on both magnesian olivines and pyroxenes. These structures, apparently of 
the orthorhombic rare-earth orthoferrite type [ Mar e zlo et al . , 1970) will be 
loosely termed perovskltes. Densities versus compositon for the perovskite 
phases and assemblages of pyroxene (ABO3) and olivine (ABO^ + AO) derived 
from l.lu's (1975a, h| X-ray parameters are also shown In Figure 7. It Is 


interesting to note that an olivine stoichiometry actually results in a lower 



12 . 


density perovsklte assemblage than a pyroxene stoichiometry. This Is a direct 
consequence of the low density o f (Mg, Fe)0 In the ABO^ ♦ AO assemblage. 

The dens It ins of such (Mg, Fe)0 recovered from high pressure lie directly 
between the densities of MgO (perlclase) and FcO (stoichiometric wilstlte) as 
can be seen In Figure 7. This Is compatible with the shock-wave data on MgO 
which shows no evidence for phase transformations to pressures over 100 GPa 
[Carter et al. , 1971; Mc Quee n and Marsh, In Birch, 1966; McQue en, 1968). 

In order to compare the static results with the shock wave data, theoretical 
metastable Hugoniots were constructed for perovsklte of appropriate composition 
(Mg. gFe, j S 10 3 ) based on Blrch-Murnaghan principal adlabats. This approach 
Is philosophically and computationally similar to that of Davies & Ga ffney 
[1973]. Necessary parameters are shown in Table 3 along with the values used 
here. Density and * -4 'e energy of the pyroxene-perovsklte transformat Ion (at Pg ) 
were derived fro: s work via Figure 7 and by using AK tr - (-PAV ♦ 1AS) (r . 

Since* the last term is not known it is assumed to be negligibly small 
(i.e., — ■ 0 for the phase boundary). Because the experimental data upper 
bound the equilibrium pressure of transition, and most likely TAS - 0, the 
resulting estimate for AK tr Is probably a strict upper bound (and rather 
high). Incidentally, though we discuss the direct transformation from pyroxene 
to perovsklte structures under shock, this obviously represents a metastable 
phase boundary (based on our current knowledge from Liu's work). Intermediate 
structures may well be involved. 

Estimates of the bulk modulus Kg and Kg' ■ / ) were based on their 

\ dl 7 P-Q 

empirical correlation with density, with particular emphasis given to data for 
perovskite-struetured compounds [D avie s, 1976; Be attie and Samara , 1971], 

The wide range In values given easily overlap estimates based on the various 
K-p relations that have been proposed [0. L. A nde rson , 1972; D. L. An derson. 

1969, 1967; D. L. An de rson a nd 0. L. Anderson, 1970; 0. L. Anderso n et al. , 1968]. 
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Similarly, GrUnei sen's (bulk or thermodynamic) y was estimated from typical 
values for various compounds, with a volume dependence that has been found 
to be not unreasonable for those cases where it has been studied. Temperature 
dependence was assumed to have a negligible effect and, again, rather wide 
bounds were given for reasonable values. Finally, a family of theoretical 
perovsklte Hugoniots was constructed for the best estimate values given In 
Table 3 and by varying each of the parameters within the stated bounds. An 
envelope containing this family of curves along with the Hugonlot based on 
the best values are shown in Figure 5. 

It Is Immediately clear that the shock-wave data Is compatible with Mg- 
pyroxene transforming to a perovskite structure at pressures above approxi- 
mately 80 GPa. This Is contrary to the conclusions of Slma kov and Trunin 
[1973). In fact, by cnsldering the most porous data (which achieve high 
densities at relatively low pressures) and *ht arguments presented above. 

It seems that perovsklte-like densities can be achieved by about 50 GPa. 

If anything, the left half of the envelope In Figure 5 agrees with the data 
as well as the best estimate, which could be consistent with an even lower 
pressure, fc * transforation to the perovskite structure. This, then, is an 

upper bound for the pressure of transformation based solely on Hugonlot data, 
yet still represents a significant overdriving compared to the static experi- 
ments [Liu, 1975b, 1976). 

Compared with the theoretical Hugoniots, the release paths still appear 
to have very steep slopes, and the most porous data cross over to yet higher 
densities. These latter data may be erroneous [Slmak ov and T runin , 1973), 
but assuming that they are not they suggest, along with the release adiabats, 
th-- possibility of a phase (or assemblage) with a zero-pressure density even 
higher than perovskite. Transformation to such a "post-perovski te" phase 
could occur at pressures as low as about 60 GPa, according to the shock data; 
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this would therefore be an upper limit for the equilibrium transition pressure. 
An initial density of 5% or so greater than for orthorhombic perovsklte could 
be consistent with the data. 

In similar fashion, a theoretical Hugoniot (Figure 6) was constructed 
for a perovskite assemblage corresponding to olivine stoichiometry (MgSi 03 
(perovsklte) + MgO (periclase)) based on the best-estimate Hugoniot in Figure 
5 the rather well determined Hugoniot of MgO [ Carter et al . . 1971; McQuee n, 
1968; McQueen and Marsh, In Birch , 1966] shown in Figure 8. Because the left 
half of the envelope in Figure 5 appeared meat consistent with the pyroxene 
points a corresponding band is shewn in Figure 6. 

In this case, perovskite -like densities are achieved in the 60-80 GPa 
range. Once more, the shock wave data are consistent with the static results 
which demonstrate olivine transferring to perovsklte at high pressures (but 
again at pressures significantly lower than under shock). On the other hand, 
the olivine points are clearly shifted to higher densities than the perovsklte- 
assemhlage Hugoniot at pressures above about 70 GPa. Because the theoretical 
curve for the perovskite assemblage cannot be shifted without causing a sig- 
nificant deviation f.om the pyroxene points, the olivine data are in strong 
disagreement with a perovskite assemblage at high pressure. 

If the olivine points did not deviate so systematically from the theoretical 
Hugoniot which is consistent with the pyroxene shock data, this divergence 
might well be ascribed to be within experimental error. However, the same 
discrepancy emerges from Figure 8, quite apart from any theoretical curves. 

The point here Is that above about 70 GPa, the olivine and pyroxene data 
overlap, whereas a perovskite model would predict olivine densities to be 
lower than pyroxene densities (i.e., shifted towards the MgO Hugoniot) as 
was mentioned above and can be inferred from Figure 7. Granted that we have 


discussed kinetic effects and a mixed phase region to the 100 GPa range, hut 
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the overlap continues to much higher pressures. Furthermore, the nature of 
the perovsklte-formlng transformations would suggest that pyroxene (which 
Involves no disproportionation) would form perovsklte at least as readily as 
olivine, whereas In the high-pressure data of Figure 8 one must either con- 
sider the olivine too dense or the pyroxene not dense enough to agree with 
a self-consistent perovsklte model. 


If*. 


CONCLUSIONS 

After reviewing experimental techniques, we p esented new Hugonlot data 
for Bamle bronzlte including simultaneously measured release paths, based 
on the mirror-buffer technique. Calculation of Hugonlot states for the 
mirror glass provides a satisfactory Internal check, with the soda-lime 
glass behaving in a fashion similar to fused quartz under shock. 

Reexamining available shock data for magnesian pyroxenes ar.d olivines 
leads to the c aclusion that they define a mixed phase (or disequilibrium) 
region to about the 100 GPa range, related to the kinetics of phase trans- 
formation in these silicates. By recognizing this point, certain discrepan- 
cies in previous interpretations of shock data can be explained. A set of 
theoretical Hugoniots for pyroxene and olivine stoichiometry, perovskit e-bear lng 
assemblages was constructed based on their properties deduced from high-pressure 
work, showing that the shock data is compatible with transformations to perovskites 
In the 45-70 GPa region (decidedly above the equilibrium pressure of trans- 
formation) . 

A perovskite model implies different densities for pyroxene and olivine 
stoichiometries, however. We note that considering an extremely simplistic 
assemblage for the lower mantle, consisting of high-pressure phases of only 
olivine and pyroxene, varying the proportions of these components results in 
a significant variation in the overall density of the assemblage. In other 
words, any compelling evidence for density variations within the lower mantle 
(based, say, on b-j ismological data) can be readily explained without resorting, 

for example, to variations in Mg:Fe ratios. This degree of freedom is 
likely to exist whatever the pertinent, multi-phase assemblage may be for 
the lower mantle. In fact, we consider our two-phase model simplistic if 
for no other reason than that we have ignored (for lack of data) the 
effects on phase relations of other likely components such '.a A1 or Ca. 
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In addition, the high-pressure shock data suggest the presence of a yet 
higher density phase (or assemblage) than perovsklte. In particular, the 
results on olivine diverge from predicted estimates and from densities con- 
sistent with pyroxene Hugo...ot points. Taken together, the data do not 
preclude transformations Involving a "post-perovsklte" phase of Mg 2 S10t, at 
pressures above 80-100 GPa, and including a reaction of the type: Mg 2 Si 2 0 6 ► 

Mg 2 Si0i, ("post-perovsklte" phase) + S10 2 (rutile, .i-Pb0 2 , or fluorite struc- 
ture: c.f., German et al. , 1974; Jamieson , this volume). 

We reiterate that the static high-pressure results (and thus our analysis) 
involve an orthorhombic (i.e., non-ideal) modification oi perovskite. Since 
this class of structures Includes a wide variety of related modifications 
[succinctly summarized by Salje , 1976], it is quite possible that the higher- 
pressure phase suggested by the shock data in fact is a perovskite such as the 
ideal cubic form. Indeed, it is quite possible that the phase quenched from 
high pressure is a modification of the actual structure attained under 
pressure. Bearing in mind the well known difficulties of recovering these 
high-pressure polymorphs on quenching [e.g., Liu and Rlngwood , 1975] and the 
important effects of environmental conditons on the perovskite modifications 
[see Sis et al. , 1973, for example] it is not unlikely that they are involved 
in higher density assemblages than have yet been described. Recently 
discussed hexagonal forms, may be an example of this [compare Burbank a nd Kvans, 

1948]. However, there are other structure-types which could also be likely, 

some of which have only begun to be discussed as possible high-pressure polymorphs 

[e.g., Moore , 1976; also, Reid and Rlngwood , 1970, and as reviewed by Rlngwood , 1975]. 

Finally, the shock data indicate very similar properties for olivine and 
pyroxene at high pressures making them both equally likely candidates for 


the lower mantle. 
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Table 1. Hugoniot Data, Soda-Lime Glass 11 * 


Shot 

No. 

Flyer 

Plate 

Projectile Velocity 
(km/sec) 

Shock Velocity 
(km/sec) 

Pressure 

(GPa) 

Density 
(gm/ cm 1 ) 

LOG 9 

Cu 

5.76 t 

0. 18 b 

8.90 

t 0.13 

91 .9 

4 . t>6 

LCC11 

Cu 

5.85 t 

0. 02 b 

8.91 

t 0.12 

93.4 

4.72 

LGC14 

Ta 

5.541 t 

o.oio c 

9.094 

♦ 0.013 

98.1 

4.756 

LtiG 1 2 

Ta 

5.66 * 

o.io c 

9.11 

♦ 0.02 

100.5 

4.85 


a) Na 2 0, 0.1; MgO, 0.02; S10 2 , 0.75; CaO, 0.13 
Initial Density, 3.49 t 0.01 gm/crn 3 

b) 1 mm-thick sample 


c) 2 mm-thick sample 


Release Adiabat 


Data. Barale Bronzite 


Table 2. 




Initial Density 

hot No. 

Flyer Plate 

(gm/cm 3 ) 

LGG6 

Cu 

1.2923 t 0.0003 

LGG8 

Cu 

3.2943 i 0.0012 

LGG14 

Ta 

3.3021 i 0.0007 


Hugonlot and 


Flyer Plate 
Velocity 
(km/ sec) 

Shock Velocity 
(km< sec) 

Pressure 

(GPa) 

5.375 a 

10.08 t 0.13 

112.6 

6.09 ♦ 0.13 

10.37 t 0.05 

132.4 

5.34 t 0.01 

10.07 t 0.03 

133.5 


Shock 
Density 
(gm/ cm') 

Buffer 
Shock 
Velocity 
(km/ sec ) 

Release 

Pressure 

(GPa) 

Release 

Density 

(gm/cm 3 / 

4.963 

9.14 ± D. 11 

80.3 

4.950 

5.261 

9.62 t 0.16 

92.3 

5.260 

5.491 

8.83 t 0.15 

73.0 

5.258 


a) Projectile velocity 


inferred from shock velocity through arrival mirrors 


Table 3. Parameters for Theoretical Hugonlot of (Mg, Fe) S10 3 Perovsklte 


Composition: (Mg_ qFe. j )S10 3 


p - 4.20 gm/cm^ 

0 

(1) 

Kg - ^260 ( + 35) GPa 

(2) 

K’ - %3.5 (±1.0) 

(2) 

AE. < ^166 KJ/mol 
tr 

(1) 

y Z 1.5 (to. 5) 



0 


>-y 0 (v/v 0 ) n 

n i 1.0 00.5) 


(1) From data of Liu (1975a, b]. 

(2) See Davies (1976), for example. 


FICURK CAPTIONS 


Figure 1. 

Figure 2. 

Figure 3. 
Figure 4. 
Figure 5. 


Figure 6. 


Figure 7. 


Flash radiographs of projectile taken at stations 390 mm (upper) 
and 33 mm (lower) from target. Nominal exposure times are M 5 nsec. 
Projectile speed is 5.5 km/sec. Shot LGG14. 

Static (a) and dynamic (b) streak camera photograph of bronxite 
Hugonlot and release adlabat experiment. Shot LGG14. 

Shock velocity versus particle velocity data for soda-lime glass. 
Shock pressure versus density for soda-lime glass. 

Compilation of shock data of Mg-rich pyroxenes. Hugoniot points 
for samples of varying porosities (synthetic, polycrystall i..e 
aggregates) and pyroxene-rich (>902) rocks are included, as are 
newly determined release paths (arrows). Theoretical Hugoniot 
(solid curve) and envelope (dashed) are discussed in the text and 
Table 3. McQueen and Marsh j "• *>66 ] are quoted in Birch (19661. 
Compilation of shock data for Mg-rich olivines including results 
on porous samples (synthetic, polycrystalline aggregates) and 
olivine-rich (>902) rocks. Theoretical Hugoniot (solid curve) 
for a Mg-endmember , high-pressure assemblage are discussed in the 
text. Dashed curve corresponds to the left side of the envelope 
in Figure 5. McQueen and Marsh (1966) are quoted in Birj’Ji (1966]. 
Summary of zero-pressure densities for high-pressure phases of 
olivines and pyroxenes as a function of composition based on 
interpretation of shock-wave data (symbols), and calculated from 
X-ray parameters for ABO3 (perovskite) and AO quenched from high- 
pressure diamond-anvil experiments (lines through points). MgO 
(perlclase) and FeO (stoichiometric wilstite) values are from: B£nard 

[1954], Katsu ra et al . [19671, Robls et al. (1966], and Rob is 
and Waldboum [1968], with error bars showing variations in measure- 
ment and estimate, while error bars on the remaining points are 
based on Liu *s [1975a, b] estimates of error in the X-ray parameters. 


Figure 7. The dashed line (ABOj + AO) corresponds to the density of a 
(Cont * d) 

perovskite assemblage with olivine stoichiometry (assuming 


Mli" p€? 

D ol perov "I)* Symbols are keyed as follows: 

A Ahre ns [1971] 

AAR Ahrens , Anderson, and Rlngwood [1969] 

AG Ahrens and Gaffney [1971] 

AS Al'tschuler a n d Sharp ldzhanov [1971a] 

I)A Davies and And erso n [1971] 

DC Davies and Gaff ney [1973] 

KP K alinin and Pan'k ov [1974] 

KHF McQ ueen, Marsh , and F ritz [1967] 

PK Pan'kov and Kali nin [1975] 

Figure 8. Representative Hugoniot data from Figures 5 and 6 compared with 

the Hugoniot for MgO [C arte r et at., 1970; McQueen , 1968; McQuee n 
and Marsh, quoted in Birch, 19661. 


1). Symbols are keyed as follows: 



